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'  Abstract 

Solutions  of  the  nonlinear  Boltzmann  equation  are  constructed  up  to  the  first  appearance  of 
shocks  in  the  corresponding  fluid  dynamics.  This  construction  assumes  the  knowledge  of  solutions 
of  the  Euler  equations  for  compressible  gas  flow.  The  Boltzmann  solution  is  found  as  a  truncated 
Hilbert  expansion  with  a  remainder,  and  the  remainder  term  solves  a  weakly  nonlinear  equation 
which  is  solved  by  iteration.  The  solutions  found  have  special  initial  values.  They  should  serve  as 
“outer  expansions”  to  which  initial  layers,  boundary  layers  and  shock  layers  can  be  matched. 


Introduction 

The  Boltzmann  equation  of  kinetic  theory  gives  a  statistical  description  of 
a  gas  of  interacting  particles.  An  important  property  of  this  equation  is  its 
asymptotic  equivalence  to  the  Euler  or  Navier-Stokes  equations  of  compressi¬ 
ble  gas  dynamics,  in  the  limit  of  small  mean  free  path.  We  propose  that,  as 
one  aspect  of  this  asymptotic  relationship,  the  question  of  existence  of 
solutions  of  the  Boltzmann  equation  can  be  reduced  to  the  existence  problem 
for  the  gas  dynamic  equations.  Although  the  latter  problem  has  received  only 
partial  solutions,  the  gas  dynamic  equations  are  much  simpler  than  the 
Boltzmann  equation  and  have  been  studied  extensively. 

This  paper  makes  a  first  step  in  that  reduction  of  the  existence  question  by 
showing  that  any  smooth  solution  of  the  Euler  equations  can  be  used  to  make 
a  corresponding  solution  of  the  Boltzmann  equation.  It  does  not  address  the 
difficulties  of  initial  values,  shocks,  and  boundaries.  The  solution  produced 
here  comes  from  a  truncated  asymptotic  expansion  and  has  special  initial 
values,  is  periodic  in  space,  and  is  valid  only  until  the  first  appearance  of 
shocks.  It  can  serve  as  an  “outer  solution"  to  which  special  solutions,  in  initial 
layers,  boundary  layers,  and  shock  layers,  can  be  matched. 

However,  the  analysis  of  these  layers  is  not  complete.  Grad  has  given  a 
formal  treatment  of  the  initial  layer  matching  problem  in  [6].  It  requires 
solving  the  spatially  homogeneous  nonlinear  Boltzmann  equation  and  showing 
that  the  solution  approaches  a  Maxwellian  asymptotically  in  time,  for  any 
initial  data.  Arkeryd  [1]  has  shown  that  this  is  true,  but  with  only  weak  L' 
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convergence,  which  is  not  strong  enough  for  the  completion  of  the  expansion. 

''I'fie  profile  of  a  weak,  steady  shock  was  analyzed  by  Nicolaenko  [11]  for  the 
Boltzmann  equation  for  hard  spheres.  The  effects  of  time  dependence  must 
also  be  included  for  the  shock  layer  analysis.  For  boundary  layer  problems, 
one  must  solye.the  steady  nonlinear  Boltzmann  equation  in  a  half -space  and 
show  approach  tljp  a  Maxwellian  at  infinity.  Guiraud  has  analyzed  a  more 
general  weakly^nonlinear  problem  in  [9]. 

fte  results  presented  here  are  valid  for  a  physically  interesting  time 
period,  until  the  first  occurrence  of  shocks,  and  concern  solutions  which  are 
far  from  spatial  equilibrium,  so  that  the  corresponding  gas  dynamics  is 
strongly  nonlinear.  Previously,  Glikson  [4]  and  Kaniel  and  Shinbrot  [10] 
showed  existence  locally  in  time.  Global  existence  of  solutions  which  start  off 
close  to  absolute  equilibrium  has  been  proved  by  Nishida  and  Imai  [12]  and 
Shizuta  and  Asano  [14],  The  asymptotic  equivalence  of  the  Boltzmann  and 
the  gas  dynamic  equations  was  demonstrated  by  Grad  [8]  and  Nishida  [13] 
for  initial  data  near  to  global  equilibrium. 

The  basis  for  the  present  work  is  the  paper  by  Caflisch  and  Papanicolaou 
[2]  on  the  Broadwell  model  of  the  Boltzmann  equation.  The  main  difficulty  in 
extending  those  results  is  the  treatment  of  the  high  velocity  tail  of  the 
distribution,  as  described  in  Section  4.  The  key  to  both  these  papers  is  the 
fact  that,  after  assuming  that  the  nonlinear  fluid  equations  can  be  solved,  the 
remaining  problem  is  only  weakly  nonlinear. 

In  Section  2,  the  existence  theorem  is  stated.  Its  proof  occupies  the  rest  of 
the  paper.  The  solution  is  found  as  a  Hilbert  expansion  with  remainder,  as 
described  in  Section  3,  and  the  equation  for  the  remainder  term  is  decom¬ 
posed  into  low  and  high  velocity  components  in  Section  4.  This  uses  a  global 
Maxwellian  distribution  wM,  as  well  as  the  local  Maxwellian  o>  associated  with 
the  Euler  equations.  Linear  and  nonlinear  collision  operators  involving  these 
two  Maxwellians  are  defined  in  Section  4.  Basic  estimates  for  these  operators 
are  presented  in  Section  5.  Then  a  linearized  version  of  the  decomposed 
remainder  equations  is  analyzed  in  Section  6.  Finally,  the  nonlinear  equations 
are  solved  by  iteration  in  Section  7. 

Use  of  the  Hilbert  expansion  in  an  existence  theorem  for  the  Boltzmann 
equation  was  suggested  by  George  C.  Papanicolaou.  I  also  want  to  thank 
Harold  Grad  for  a  number  of  discussions.  Most  of  this  work  was  completed  at 
the  Mathematics  Research  Center  and  the  Courant  Institute  of  Mathematical 
Sciences,  whose  support  I  am  happy  to  acknowledge. 

2.  The  Existence  Theorem 


The  Boltzmann  equation  is 


(s+«'v)f4cw;fi 


f.-rwr---  v'nv)'*.' 
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in  which  V  =  3/ax.  The  function  F  =  F(x,  t,  £)  is  the  density  of  particles  of 
velocity  | €  J?3,  at  position  xe[0, 1]  and  time  re[0,  °°).  We  look  for  solutions 
which  are  spatially  periodic  in  [0, 1],  The  left-hand  side  of  (2.1)  represents 
streaming,  while  Q  is  a  nonlinear  integral  operator  representing  collisions. 
Following  Grad  [6],  we  consider  cut-off,  hard  potentials  only.  The  parameter  e 
is  proportional  to  the  mean  free  time  and  is  assumed  to  be  small. 

There  are  special  distribution  functions  to,  called  Maxwellians  and  given  by 

(2.2)  <•>(€)  =  exp  {-(£- u)2/2  Tj , 

which  are  in  equilibrium  with  the  collision  process,  i.e.. 


(2.3)  Q(&>,  to)  =  0  . 

If  p,  ii,  T  are  constant  in  x  and  t  (as  well  as  in  |),  a>  is  called  a  global 
Maxwellian :  if  they  depend  on  x  and  t,  it  is  a  local  Maxwellian.  The  constants 
p,  u,  T  are  the  macroscopic  density,  velocity  and  temperature,  respectively. 
More  explicit  descriptions  as  well  as  derivations  and  basic  properties  of  the 
Boltzmann  equation  can  be  found  in  [5]. 

The  fluid  dynamic  description  of  a  gas  is  given  by  the  Euler  equations 

|-p  +  Vpu,  =  0, 

(24)  ^pB  +  V(pu1n)  +  Vp  =  0, 

~  p(e  +2M2)  + V(pu,(e  +|h2))  + V(puj)  =  0 , 
p  =  pRT  =  | pe . 

These  are  the  equations  of  conservation  of  mass,  momentum,  and  energy,  and 
the  equation  of  state. 

The  following  theorem  shows  that  there  are  solutions  of  the  Boltzmann 
equation  which  are  nearly  local  Maxwellian,  in  which  the  macroscopic 
variables  p,  a,  T  evolve  according  to  the  Euler  equations. 


Theorem.  Let  (p(x,  t),  a(x,  t),  T(x,  t))  be  a  smooth ,  spatially  periodic 
solution  of  the  Euler  equations  (2.4)  for  t  e  [0,  t],  x  e  [0, 1].  Construct  the  local 
Maxwellian  w(x,  f,  £)  from  p,  u,  T  as  in  (2.2).  There  is  a  positive  e0  such  that , 
for  each  0  <  e  S  e„,  a  smooth  solution  F *  of  the  Boltzmann  equation  (2.1)  exists 
for  t  e  [0,  t],  with 

t  DTIC  TAB 

F*  e  L~([0,  T];  H'(x,  {))  n  C([ 0,  TJ;  L2) , 


(2.5) 


-F*eL~([0,  r];  L2) . 
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Moreover, 

(2.6)  UP  -  (o\\  S  ce , 

where  the  norm  is  in  any  of  the  spaces  in  (2.5)  and  c  is  independent  of  e. 

L“([0,  t];  L2)  means  the  Banach  space  of  bounded  measurable  functions 
from  [0,  t]  to  L2([0, 1],  R3),  etc. 

Remarks.  (1)  This  theorem  assumes  a  solution  of  the  nonlinear  Euler 
equations  is  in  hand.  Although  these  equations  are  much  simpler  than  the 
Boltzmann  equation,  the  existence  question  has  been  only  partially  answered. 

The  smoothness  of  the  Euler  equations  is  expected  to  last  until  the  first 
appearance  of  a  shock,  which  we  would  take  as  the  time  t.  We  have  not 
optimized  the  amount  of  smoothness  needed,  but  certainly  it  is  enough  for  p, 
u,  T  to  be  in  H9. 

(2) .  The  initial  values  of  P  are  essentially  those  of  the  Hilbert  expansion 
(see  Section  3).  This  is  just  what  is  needed  to  match  to  an  initial  layer. 

(3)  The  one-dimensional  spatially  periodic  problem  is  handled  for  simplic¬ 
ity.  We  could  just  as  well  do  the  multi-dimensional  infinite  space  problem, 
which  is  needed  for  the  matching  to  boundary  layers. 

(4)  The  Navier-Stokes  equations  could  be  used  instead  of  the  Euler 
equations.  The  viscosity  is  multiplied  by  e  and  uniform  smoothness  as  well  as 
nice  limiting  behavior  is  required  of  the  solutions.  The  result  is  an  approxima¬ 
tion  with  error  size  e2. 

3.  The  Hilbert  Expansion  with  Remainder 

The  solution  P  is  found  as  a  sum 

(3.1)  p=  i  e"Fn  +  e3PR, 

where  F0,  •  •  • ,  F6  are  independent  of  e.  They  are  the  first  6  terms  of  the 
Hilbert  expansion,  which  solve  the  equations 


(3.2) 

0=Q(Fo,Fo), 

(3.3) 

(^+«i  v)fo=20(F0,F,), 

(3.4) 

(£  +  *iV)f,  =  2 Q(F0,  F2)+  Q(F, ,  Ft) , 

(3.5) 

(£+*,v)f5  =  2Q(F0,F6)+  I  Q(Fj, Fj) 

7  i+j= 6 

13»S5 
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From  equation  (3.2)  we  infer  that  F0  =  o>,  a  local  Maxwellian.  The  remaining 
equations  (3.3M3.5)  involve  the  linear  operator  if  =  -2Q(F0,  * ),  an  integral 
operator  over  the  velocity  space  (|).  This  Fredholm  operator  can  be  inverted 
after  checking  that  the  inhomogeneity  is  perpendicular  to  the  null  space 
{</>o ,  •  •  • ,  4>s)  of  the  adjoint  operator  i£t.  The  appropriate  inner  product  is 

</.g>=(  /(6)g(«)d|, 

Jr3 


and  the  are  given  by 


(3.6) 


4>o  =  ~ , 
P 

pT 


so  that  <£,)  =  5i(.  The  null  space  of  if  is  {o*£0,  •  •  ■ ,  ax£s}. 
Equations  (3.3)  and  (3.4)  become 


(3.7) 

(3.8,  + 

(3.9) 

(3.10)  F2=-if-,((|  +  ^v)F1-Q(F1,F)))  +  «l>2, 


in  which  if<t>,  =  if<62  =  0.  Equation  (3.7)  gives  exactly  the  nonliner  Euler 
equations  (2.4)  for  the  macroscopic  density,  velocity,  and  temperature  p,  n,  T 
corresponding  to  <v  =  F0.  From  equation  (3.8),  we  see  that  F,  = 
(pl4>0+u’<hi  +  T'd>4)ci>  +  '|'1  in  which  if ’¥,  = -(3/dt  +  |,V)F0  and  'P,€N(ify\ 
The  coefficients  pl,  n\  T1  are  the  macroscopic  density,  velocity  and  tempera¬ 
ture  corresponding  to  F] .  According  to  (3.9)  they  satisfy  inhomogeneous 
Euler  equations  linearized  about  p,  n,  T  and  with  inhomogeneity  given  as  an 
operator  on  p,  a,  T.  The  terms  F2,  •  ■  • ,  F6  are  found  similarly.  However  since 
the  expansion  is  truncated  we  can  take  p6  =  u6=  T6  =  0. 

A  careful  treatment  of  the  Hilbert  expansion  is  found  in  [5]  and  [7].  Only 
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several  facts  are  needed  here.  We  are  starting  with  a  smooth  solution  (p,«,  T) 
of  the  nonlinear  Euler  equations,  from  which  we  construct  F0  =  o>  which 
solves  (3.7).  The  remaining  Hilbert  expansion  equations  are  linear  and  have 
solutions  which  are  smooth  in  (x,  t)  and  decay  in  fc  Consider  F, .  Grad  has 
shown  in  [7]  that  preserves  decay  in  £,  so  that  ~ 

|||3  exp{-(|-n)2/2T}.  Since  the  inversion  is  local  in  (x,  1),  'P,  is  smooth  in 
(x,  l).  The  remaining  terms  in  F,  obviously  decay  like  |||2  exp{-(|-«)2/2T},  and 
they  are  smooth  in  (x,  r)  since  coefficients  p1,  u1,  T1  solve  linear  equations 
with  forcing  terms  coming  from  the  smooth  functions  p,  a,  T.  The  following 
proposition  summarizes  these  facts. 

Proposition  3.1.  Let  (p,  o,  T)  be  a  smooth  solution  of  the  Euler  equations 
(2.4),  and  form  the  Maxwellian  F0  =  to  as  in  (2.2).  Then  the  terms  F, ,  •  •  • ,  F6 
of  the  Hilbert  expansion  are  smooth  in  (x,  f)  and  have  decay  given  by 

(3.11)  Fj  (x,  t,  |)  =  c  |||3‘<«j(x,  t,  |) , 
where  c  is  a  constant  independent  of  |,  x,  t. 

Next  we  find  an  equation  for  the  remainder  Fr,  by  putting  the  expansion 
(3.1)  into  equation  (2.1)  and  subtracting  the  Hilbert  expansion  equations 
(3.2M3.5).  After  dividing  by  e3  and  regrouping,  the  equation  becomes 

(3. 1 2)  1 ' v)fr  =  ~  Q(a>,  Fr )  +  2Q(F,  +  eF2  +  e  2F3 ,  F„ ) 

+  e20(Fr,  Fr)  +  e2A , 


where 


(3.13) 


fs+  2*  e 

iSS.jSS 


i+j— 6. 


0(F, ,  Fj) . 


For  initial  values  we  take 

(3.14)  FR(t  =  0)  =  0. 

The  superscript  e  on  FR  has  been  dropped.  Once  the  remainder  equation 
(3.12),  (3.14)  is  solved,  we  have  the  desired  solution  of  the  Boltzmann 
equation.  By  writing  the  solution  as  in  (3.1)  we  have  reduced  the  Boltzmann 
equation  to  an  equation  in  which  the  nonlinearity  and  inhomogeneity  are 
both  small. 
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4.  Decomposition  of  FK 

In  the  past  the  linearized  Boltzmann  equation  has  been  analyzed  by 
symmetrizing  the  operator  i?=  -2Q(w,  •)  to  get  L=  -  2a>1/2Q(a>,  to' 1/2  • ). 
Since  L  is  non-negative  the  (l/e)L  term  does  not  cause  growth.  This  is  fine  if 
a)  is  constant  in  (x,  f),  but  otherwise  the  symmetrization  procedure  results  in  a 
new  term  {u>~V2(dldt  +  (idldx)(<oll2)}f  which  is  like  |£| af  and  has  an  uncontrol¬ 
led  sign.  Thus  at  large  velocities,  the  distribution  function  may  be  growing 
rapidly  due  to  streaming.  Similar  problems  occur  in  the  study  of  soft 
potentials  (cf.  [3]),  shocks  (cf.  [11]),  and  boundaries. 

To  remedy  this  difficulty,  we  decompose  FR  into  essentially  low  velocity 
and  high  velocity  parts  using  two  different  Maxwellians — to,  as  above  and 


(4.1)  <*>M  =  —  *  --7-^ exp{-g2/2Tm} , 

(Z7TlM/ 

in  which  TM  is  constant,  and  TM  >max1I  T(x,  t)  so  that  coMSca>  for  all 
(x,  t,  |).  We  shall  employ  various  operators  which  are  linearized  about  these 
Maxwellians. 

Notation. 

(4.2)  L/=  -2<a-mQ(<o,o>mf)  =  (v  +  K)f, 

where  v  is  a  function  of  |  and  K  is  a  compact  integral  operator  over  |, 

(4.3)  !>(!)=  j<o.dn, 

(4.4)  g)  =  a rmQiwlnf,  wmg)  . 

The  definitions  above  are  those  of  Grad  [7],  [8]  (the  cross  section  di 1  is 
defined  there).  The  following  are  new,  but  are  analogous  to  Grad’s  defini¬ 
tions: 

(4.5)  Lm/=  -  2wm,/2Q(w,  =  ( v  +  KM)f , 

(4.6)  L\f-  2wmI/zQ(F1  +  eF2  +  e2F3,  wtff) , 

(4.7)  g)  =  <otfg) , 

(4.8)  vM(|)=|«jMdn. 
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We  also  define 


(4.9) 

f!  for|g|^0 
10  for  |€|>^0 

(4.10) 

X  =  1-X> 

(4.11) 

a  =  , 

(4.12) 

ft  =3W*l(^  +  ^iVjo) 

(4.13) 

a  =  a>MU2A. 

Note  that  a  is  exponentially  decaying  since  TM  >  T  for  all  x.  t. 

Now  decompose  FR  as 

(4.14)  FR  =  o>1/2g  +  <i>£i2 h . 

Essentially,  <oxn g  is  the  low  velocity  part  and  <o}J,2h  the  high  velocity  part;  the 
precise  definitions  are  that  g  and  h  solve 

(4.15)  (£  +  €tv)g=-jLg-±x<r-'KMh, 

(4.16)  (~  +  tiV)h=~  nag  — (r  +  xK^h+Ltiag  +  h) 

+  e2vMrM(ag  +  h,ag  +  h)  +  e2a, 

(4.17)  g(t  =  0)  =  h(t  =  0)  =  0  . 

If  (4.15)  is  multiplied  by  o>1/2  and  (4.16)  by  wjf,  and  the  two  are  added,  the 
result  is  just  equation  (3.12)  for  FR  defined  by  (4.14).  After  solving  for  g  and 
h,  the  solution  F  of  the  Boltzmann  equation  (2.1)  will  be  complete. 

5.  Basic  Estimates 

Before  analyzing  equations  (4.15)  and  (4.16),  basic  estimates  are  needed 
for  the  operators  defined  in  the  last  section.  First  we  define  some  norms. 

Definitions. 

(5.1)  ll/ilr  =  sup  sup (l  +  |£|)r(  [/(*,!,  t)2dx)  , 

OStST  t  W  / 

(5.2)  \\fl,=  t  IIV7IL 

n-0 
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(5.3) 

111/(0111=  (J/(M,  02d|dx)  , 

(5.4) 

111/(0111,  s  t  |||V»/(0|||. 

n  “0 

In  the  estimates  of  this  paper  we  use  c  as  a  generic  constant.  Any  other 
constant,  e.g.,  c2+ 1,  will  be  replaced  by  c.  Often  all  constants  will  be 
omitted. 

Basic  Estimates. 

Grad  [7],  [8]  proved  the  following  estimates: 

(5.5) 

nmsii/iu, , 

(5.6) 

\\Kf\\o=  sup  Hl/lll, 

OSISt 

(5.7) 

IIIk/I||s|||/IHs||/||2, 

(5.8) 

a/./)so, 

(5.9) 

.  Iir(/.g)lls  ll/ll,  ,||glL 

rSl  , 

(5.10) 

cS„(|)S|||. 

There  are  analogous  estimates  for  the  operators  KM  and  T M : 

(5.11) 

IIKm/IIsII/IU, 

(5.12) 

Hkm/11oS  sup  m/lll , 

OStST 

(5.13) 

III^m/HI = lll/lll = ||/||2 , 

(5.14) 

lirM(/,g)is||/IUIg||„ 

rS  1 . 

We  also  need  the  following  new  estimates: 
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(5.18)  ||i(VK)/|s||/|U2, 

(5.19)  km§ci/, 

(5.20)  Vv  =  cv. 

The  integral  operator  KM  is  written  explicitly,  as  in  the  analysis  of  Grad 
[7],  in  the  form 

«*/(©=[ 

Jr' 

(5.21)  (,(€,t1)  =  27r-^5a>ii2(n)f  B(0,v)d», 

wM(l)  Jr’ 

eXp{"^^2+W)2}0(V>W)dW> 

where  v  =  n  - 1,  £  =  |(|  +*l)  =  Ci  +  £2,  with  n,  the  component  of  a  parallel  to  v, 
5,  the  component  of  £  parallel  to  v.  The  estimates' (5.11),  (5.12),  (5.13)  are 
proved  using  the  bound 

(5.22)  |  S  c|  exp  { -  c(v2  +  {*)} . 

The  estimate  (5.14)  is  exactly  (5.9)  with  co  replaced  by  com-  The  ine¬ 
qualities  for  o,  I'm,  Vo  are  derived  directly  from  the  definitions  (4.3),  (4.8). 
Bounds  on  L,,  VLt,  VLM,  VK  are  proved  using  (5.9)  or  (5.14),  (5.19),  and 

(5.20),  and  the  smoothness  and  decay  of  the  F(. 

The  next  lemma  uses  the  basic  estimates,  the  exponential  decay  of  a,  and 
the  smoothness  and  decay  of  the  Ff. 

Lemma  5.1.  The  following  estimates  hold: 

sup(l  +  |5|)V(|)Sc, 

( 


(5.23) 
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£ 

F 


(5.24) 
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ll«Xllr.«  S  C  ||g|lo.,  , 


(5.25)  l^^r^og.  +  fc,,  ag2  +  h2) |  S (llgillo..  +  l|hiD(l|g2|lo.,  +  HMD . 

rSl,  sSl, 


(5.26) 


1 

—  a  Sc. 
»»  II,, 


In  Section  6,  the  following  modified  Gronwall  inequality  is  used. 
Lemma  5.2.  Let  u(t)  5  0,  a  SO,  0SO,  and  suppose  that 


(5.27) 
Then 

(5.28) 


~u(t)  =  a  sup  u(s)  +  £(t). 
at  osisi 


u(t)S  h(0)c“‘  + 


1“'- 


>|3(s)  ds. 


6.  The  Linearized  Equation 

The  linearized  version  of  equations  (4.15)  and  (4.16)  is: 

(6.D  (^  +  |1v)g=-iLg-^a-1KMh, 

(6-2)  (ar  +  ^*V),I=  -fitr8"(*'+^KM)h+jLi(o-g  +  h)  +  e2I>, 

(6.3)  g(t  =  0)=  h(f  =  0)  =  0  . 

The  estimates  in  the  next  lemma  are  the  key  steps  in  this  paper. 
Lemma  6.1.  Let  g,  h  solve  equations  (6.1),  (6.2),  (6.3).  Then 

(6-4)  llg|lr.1Se1/4||f>/i'||,+I,t , 

(6-5)  Mr.^^llb/HL, 

for  3,  sS  1. 
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Proof:  (A)  Estimates  on  h.  The  integral  version  of  (6.2)  is 

(6.6)  h(x,  r’€)  =  |  exp{-(t-s)r'(t)/£}|-/iag-jx/CMh  +  L,(o-g  +  )j)+e2fe 
x(x-(t-s)^,,  £,  s)  ds. 

Using  Grad’s  basic  lemma,  [8],  p.  164,  we  estimate 


\Mr  =  e  II-  jMTgjl  +  II-  xKMh\\  +  e  ||i  L,(erg  +  h) ||  +  e3  ||-  b 


(6.7) 


=  «  llg|lo+  —  Ml,  +  e(lig|lo  +  Ml, )  +  £3  \\blv\\r , 

c0 


making  use  of  (5.10),  (5.11),  (5.15),  (5.24),  (4.10).  By  choosing  £0  large 
enough,  it  follows  that 


(6.8)  MrSe||gi|0+e3||&/HI,. 

We  differentiate  (6.2)  to  get  an  equation  for  Vh : 

(6-9)  +  -(M<rVg  +  (V(x<j)g)-^(i/  +  ^6.'M)Vh -^-(Vi>  +  *Km)Ji 

+  L^o-Vg  +  Vh  +  (Vcr)g)  +  (VL,)(trg  +  h)+  e2Vb , 

and  continue  as  before  to  obtain 

(6-10)  ||Vh||rSe||g||0,1  +  e3||b/4+£3jivb|  . 

Combined  with  (6.8),  this  says  that 
(6.11)  IIMUiSe  ||g||0.1+e3  ||6/l4i 

(B)  Estimates  on  g.  The  integral  version  of  (6.1)  is 
g(x,t,|) 

=  J  exp  { -  (f  -  s)i/(|)/e}|  -i  Kg  -  j  x<r~'  KMh|(x  -  (f  -  s)£, ,  (,  s)  ds . 
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Making  use  of  Grad’s  lemma  we  find 

(6-13)  ||g|lrS||Kg||r+|Uo--‘KMh||r.  ' 

But 

||<r_1xXMlj||r  £exp{c&>} 

(6.14)  =  exp  {c^oKe  l|gllo+e3||f)/Hlr-i) 

=  ein  ||g||o  +  es/2 IlftMl-i . 

using  (5.11)  and  (6.8)  and  choosing  e  small  enough  so  that  e  exp{c£o}<£I/2 
The  first  term  in  (6.13)  is  estimated  by  (5.5)  or  (5.6),  with  the  result  that 

(6.15)  Ilgllr^llglLi  +eS/2l|6/HL-i .  rgl, 

(6.16)  ||g||0S  sup  ||ig|||  +  e5/2||f>Mlo. 

OS!  St 


Employing  (6.15)  for  r  and  then  for  (r-1),  recursively,  and  finally  using 

(6.16) ,  yields 

(6.17)  l|g||r  =  sup  lllglll-t-  e5/2  llfo/Hlr-1  - 

OSlST 

To  estimate  |||g||],  we  multiply  equation  (6.1)  by  g  and  integrate  over  x  and 
£  We  use  the  spatial  periodicity  and  the  non-negativity  of  L  to  drop  the  d/dx 
and  the  L  terms,  the  result  is 

(6.18)  “ll|g(»)lll2s1llk-,^Mh||Hllg|||. 

2  of  s 

Thus 

~  lllg(0l!l  =  exp  {c£o)  -  ||h||2 

at  S 

(6-19)  S  exp  {c|oKi|g||o  +  e2  \\bl  Hlz) 

=  exp{c^}(  sup  |||g(s)|||  +  e2||6/H|2). 

Since  this  is  true  for  any  t,  t St,  the  modified  Gronwall  inequality  is 
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applicable  and  implies  that 


(6.20) 


|||g(r)|||  S  e2  ||b/i4  exp  {t  exp  {cfjj}} 

Se3/2||b/HK 


by  choosing  e  small  enough  so  that  ve exp {r exp {c£o}}  =  1.  Substitution  into 
(6.17),  gives  the  result 


(6.21) 


||g||,sie3'2IMU, 


rS3 . 


(C)  Estimates  on  Vg.  The  equation  for  Vg  is 

(^+SiV)vg  =  -~(v+  K)Vg--(Vv  +  VK)g 
\dt  /  e  e 

(6.22)  -  -  x(V(<t-,)A:m  +  a'  l(VKM  ))h--x<r~'  KM(Vh ) , 

£  £ 

and  proceeding  as  before  we  find 


(6.23)  l|Vg||r£  sup  |||Vg||[+e3/2||f>M|r+1 ,  . 

OStST 


The  L2  estimate  on  Vg  is  found  by  multiplying  equation  (6.22)  by  Vg  and 
integrating  over  x  and  As  before,  the  x  derivative  and  L,  g  terms  can  be 
omitted.  Hence 

“ll|Vg|||2S--((V^+VK)g,Vg) 

Z  at  £ 

(6-24)  -  -  U(V(a-,)KM  +  a-\VKM))h,  Vg) 

e 

-i(xa-‘KM(Vh),Vg). 

£ 

We  use  the  Schwartz  inequality  to  estimate 

£  IllVgtll  S -  |||(Vp  +  VK)g|||  +  -  ||MV(«7-')Km  +  <r-\VKM))h\\\ 
dt  e  e 

+  -  lllx^-1  Vh|||  s-  ||g||4  +  -  exp  {c|g}  HMki 

£  £  £ 

=  e 1,2  llb/HIa  +  exp  {cfgKKgllo.1  +  e2  ||6/i4.,) 

Sexp{c|l}  sup  i||Vg|||  +  e1/2||b/F||3il . 

Qatar 


(6.25) 
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From  the  Gronwall  inequality  we  obtain 

(6.26)  |||Vg|||(r)Se 1/2 \\b/ Hki  exp{exp(c£o}}Se,/4||b/i4.i , 

by  choosing  e  sufficiently  small.  Combining  this  with  (6.23),  (6.21),  (6.8),  and 
(6.11),  we  get 

l8lL.1S«M4IWHU,,  '*2, 

IMr.,Ses,4||ty«4.i.  r  =  3 . 

The  same  estimates  are  true  with  s  >  1 ;  thus  the  proof  of  Lemma  6. 1  is 
complete. 


7.  Solution  of  the  Nonlinear  Equations 

The  nonlinear  equations  (4.15),  (4.16),  (4.17)  are  solved  by  the  implicit 
function  theorem  or  direct  iteration  on  the  linear  equations.  Let  g„+J,  h„+1  be 
the  solution  of  equations  (6.1),  (6.2),  (6.3),  with 

b  =  vMYM(crgn  +  hH<  (Tg„+hn)  +  a, 


and  g0  =  h0  =  0.  Then 

(7.D  II- fill  s»gJ&i+|JUEl+cf 

II v  Hr.  1 

and  from  Lemma  6.1,  we  get 

Il&i  +  lllr.l  =  eI/4(l!&t!lo,l  +  llbnlloi  +c)  , 

'  lk+1|i,1Sew(||g„|g.1+lkll?.1+c). 

By  iteration  we  find 

(7.3)  »&.+iH.iSc,  H^.llr.tSc, 

in  which  c  is  independent  of  n.  The  differences  gn+]-g„,  hn+,  -  hn  solve 
equations  (6.1),  (6.2),  (6.3)  with 

b=  r-rnTM^gn  +  fin+ergn,^^-,,  Og*  +  K  ~  0%,-,  ~  K-\)  > 

|i  b|  —  c(||gn  -  gn—  ljlo.l  +l|bn  -  Hh-iIL.i)  . 


(7.4) 
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so  that 

_  5)  llgn  +  1  -  gjr.l  S  el/4C(||g„  -  8*.  —  1  Ho.  1  +WK  ~  hn-lllr.l)  , 

ll^n+1  -  M,.l  ^  £  ‘/4c(Hg„  ~  +  IK  “  - 1  lir.  1 )  • 


From  this  it  follows  that  g„  — »  g,  h„  —»  h  in  ||  •  ||M .  By  standard  arguments  one 
then  finds  that  g,  h  are  solutions  of  (4.15)  and  (4.16).  This  completes  the 
proof  of  Theorem  2.1. 
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